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^ ■ ABSTRACT: The exact master equation for a harmonic oscillator coupled to a 

heat bath, derived recently by Hu, Paz and Zhang, is simplified by taking the weak- 
, coupling, late-time limit. The unique time- independent solution to this simplified 

^ . master equation is the canonical ensemble at the temperature of the bath. The 



frequency of the oscillator is effectively lowered by the interaction with the bath. 
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INTRODUCTION 

The evolution of the density matrix of a quantum harmonic oscillator linearly 
coupled to a heat bath is a fundamental problem, but it has only recently been 
solved exactly [1]. Physical quantum mechanical systems have often been idealized 
as isolated, with the sole effect of the environment being the maintenance of a finite 
temperature. The non-trivial role of the environment in decohering the excitations 
of a weakly-coupled system is currently still being explored; as a complement to this 
research, this letter examines the ground state of a harmonic oscillator weakly coupled 
to a heat bath. 

If a system is in equilibrium with a thermal environment, to which it is weakly 
coupled, it has long been assumed that the reduced density matrix for the system 
is given by the canonical ensemble at the environmental temperature. Using the 
density matrix evolution equation ('master equation') derived in Reference [1], this 
letter confirms this assumption in the case of the harmonic oscillator. The discussion 
assumes familiarity with the results of Reference [1] (hereafter denoted HPZ), which 
will not be derived here. 



THE HPZ MASTER EQUATION AT 
WEAK COUPLING AND LATE TIME 

In HPZ, the following master equation is derived, for the time evolution of the 
density matrix p{x, x') of a simple harmonic oscillator, with the position variable x 
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coupled linearly to a heat bath: 
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" dx dx' 
Here h and the oscillator mass have been set equal to 1. The time-dependent real 

co-efficients defined in HPZ have been renamed A, B, C, and D. The terms in (1) 

proportional to A and D are responsible for diffusive effects in the evolution of p. 

B could be considered a time- dependent addition to the effective frequency of the 

oscillator, while C has the effect of a time-dependent dissipation constant. 

Equation (1) is derived using an initial quantum state that is a direct product of 
the initial oscillator and (thermal) environment states; the authors of HPZ suggest 
that some features of their master equation may be artifacts of this actually rather 
implausible initial condition. In order to avoid the spurious time dependences intro- 
duced by the artificial assumption that the oscillator and heat bath are uncorrelated 
at time t = 0, the co-efficients will be replaced by their asymptotic forms at late 
times: A{t) -^ A^ = A{oo), etc. Setting the RHS of (1) equal to zero, one obtains 
an equation for the time-independent state into which the open harmonic oscillator 
might be expected to settle down at late time. Equivalently, this equation can be 
considered to describe the ground state of a simple harmonic oscillator in the presence 
of a thermal environment: 
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The time-independent co-efficients have simple forms in the weak-couphng hmit, 
where only terms of up to second order in the bath-oscillator coupling constant are 
retained[HPZ, eq. 2.46]. In the notation of the present letter, 

Aoo =9 jz- lim ds rfu; /(a;) coth — cosct;ssinr2oS 
iZo t^oo Jq Jq 2 
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B^ = — g lim I ds I du I{ijj)sinuJSCOsVtQS 

t^oo Jq Jq 

Coo =9 TT- liiii ds dio I{uj) sinus siniloS 
"0 i^oo Jq Jq 

rt rT g^ 

Doo =9 hm ds dujI{u)coth — cosajscosl^o-s • 
*^oo Jo JO 2 

I{uj) is the spectral density of the environmental heat bath, and 9 is the bath-oscillator 
coupling constant. It is assumed that I{uj) is cut off at some high frequency F and 
vanishes at least linearly at zero frequency, so that, as long as the limit t -^ 00 is 
taken last, the order of the s and u integrals is arbitrary. The temperature of the 
bath is kT = P^ , where k is the Boltzman constant. 

In solving (2) to ffist order in g , the co-efficient -Boo will be split into two terms: 

Bo, = ilo tanh f^Aoo + f^o^^ • (4) 

The first term is chosen to provide a cancellation shown below, while the second term 
is a renormalizing correction to the effective frequency of the oscillator, proportional to 
g . {Sn will be determined in the next section of this letter.) Define the renormalized 
frequency to be fi = fio + SQ. 

A standard representation of the delta function allows one to write 



One can therefore re-write (2) in the final form 
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keeping only terms of up to first order in g . This is the weak-coupling, infinite-time 
limit of the exact time- independent master equation of HPZ. 



Define the thermal density matrix 



pp{x,x') = (1 - e-^") Y, e-^P^^^{x)^^{x') 

n=0 



(7) 



where ipn{x) is the wave function for the nth excited state of the harmonic oscillator 
with frequency Vt. It will now be verified that p = pp is a. solution to (6). The 
verification is straightforward, and proceeds line by line. 

The first line of (6) is simply [Hx — Hx']pp, for H the harmonic oscillator Hamil- 
tonian in the position representation, and so equals zero. 



The second line is proportional to 
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Here a = —X={ilx + -§^), a' = '^^{^^ ~ ^)' ^^^ primes imply that x is replaced by 



/2n^"~ ' ^^' " ^20. 

x\ and ^„ satisfies 

a'^„ = Vri + 1 ^„+i 
^_i = . 

Similarly, the third line of (6) is proportional to 
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The canonical ensemble at temperature kT = /3^ is therefore a solution to 
the late-time, weak coupling, time-independent master equation. This equation is 
a hyperbolic partial differential equation in the two independent variables x and x' . 
Normalizability of p requires that p and its derivatives decay to zero for large |x| or 
\x'\. Normalizability therefore imposes Cauchy initial, final, and boundary conditions 
on the master equation, which in general over-determine the solution[2]. The solution 
which has been found is therefore unique. (The terms "initial" and "final" are used 
by analogy, since either x or x' may be considered to play the role of the timelike 
variable usually involved in hyperbolic equations.) 
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DETERMINATION OF SQ 



From equations (3) and (4) one can determine the frequency correction 6Q. Com- 
paring the definitions of Coo and D^o, one sees that 

(^f] = -^lim [ duI{oo)coth-^R{u,n,t) , (10) 

\l i— >oo Jq 2 

where terms of order g are neglected, and R is defined by 

duj r^ BVt r^ 

i?(cj, fi, t) = tanh — / ds sina;scosr2s + tanh / ds cos uJssinVts . (11) 

R can be evaluated by expanding the hyperbolic tangents in Taylor series, and noting 
that, for integer n > 0, 
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One then has 
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where A„ are the co-efficients in the Taylor series for the hyperbolic tangent. 

The oscillating functions of t in the last two lines of (12) may in fact be ignored 
in the limit t -^ oo: since I{uj) vanishes (at least) linearly at a; = 0, /(a;)coth^ 
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may be expanded in a series of non-negative powers of w, and explicit calculation will 
show that 

lim / dxx'^e'*'' = (14) 

for m > 0. 

Ignoring the t-dependent part of Rn{uj,Q,t), one therefore has, to leading order 



in g^, 

'0 
In the high temperature limit, this becomes 
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while in the low temperature limit it approaches 
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CONCLUSION 

In recent years it has been argued that open quantum systems are generically de- 
cohered by the environment to which they are coupled, in such a way that the state 
of the system is rapidly driven towards a mixture of eigenstates of the interaction 
IIamiltonian[l,3,4,5]. This phenomenon is considered to occur on a short time scale; 
open quantum systems in equilibrium after a long period of time have received com- 
paratively little attention. It is nevertheless of fundamental importance to confirm 
that the equilibrium state of a system whose position variable is weakly coupled to a 
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heat bath is indeed the canonical ensemble, and not a mixture of position eigenstates. 
This is true at all temperatures, and is independent of the spectral density of the heat 
bath, to leading order in the coupling g . 

At higher order in g\ this result probably no longer holds. Heuristic arguments 
justifying the canonical ensemble typically assume weak coupling between the system 
and its environment. Presumably the equilibrium solution for p depends at higher 
orders in g on the specific form of the spectral density I{u>). 

The infinite-time limit of the master equation derived using uncorrelated initial 
states at t = may be conjectured to be equivalent to the master equation one 
would obtain at finite times using more realistic initial conditions. The canonical 
ensemble is an explicit example of a state which does not (at leading order in the 
coupling constant) suffer wave-function collapse onto the pointer basis determined by 
the coupling to the environment. This supports the suggestions by several previous 
researchers that the role of the initial conditions in decoherence and environmentally- 
induced superselection needs further investigation. 
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